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SUMMARY 


A model of the spanwise variation of deviation for stator 
blades is presented. Deviation is defined as the difference 
between the passage mean flow angle and the metal angle at the 
outlet of a blade element of an axial compressor stage. The 
variation of deviation is taken as the difference above or below 
that predicted by blade element, i.e., two-dimensional, theory at 
any spanwise location. 

The variation of deviation is dependent upon the blade 
camber, solidity and inlet boundary layer thickness at the hub or 
tip end-wall, and the blade channel aspect ratio. If these 
parameters are known or can be calculated, the model provides a 
reasonable approximation of the spanwise variation of deviation 
for most compressor middle stage stators operating at subsonic 
inlet Mach- .numbers.- . 


INTRODUCTION 


The development of a spanwise deviation model for stator. 1 ?., 
was done before the rotor model because the effects of rotation 
and end-wall gap could be ignored. These factors lead to a 
significantly different variation of deviation between rotors 
and stators, as will be shown. 
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A model for the spanwise distribution of deviation can be 
split into two parts: the prediction of (1) primary or two- 

dimensional blade element deviation, and (2) secondary flow or 
end-wall induced deviation. This paper considers alternative (2). 

NASA data [1]* from twelve (12) subsonic, middle stage 
rotor/stator combinations have been considered to form a model of 
stator spanwise secondary flow induced deviation. All data 
considered was from operation at the design tip speed of 244 
meters per second and maximum efficiency. A summary plot of 

•j* 

this data is shown in Figure 1. There, the measured deviation, 6 , 
minus a calculated two-dimensional equivalent deviation, 6 eg, 

(2] is plotted versus percent span. If the equivalent deviation 
was correct and the flow near mid-span was approximately axr- 
symmetrical (i.e., two-dimensional in the blade element plane), 
the mid-span values of <5 ~ <5 eq shown on Figure 1 would be zero. 

That it' is' not implies" a lack" of core flow, an error in the 
calculated blade element performance, or both. If it is assumed 
that the calculated blade element deviation is mainly at fault, 
then a replot of the data on Figure 1 with the mid-span value of 
<5 - <5eq set at zero would give the spanwise variation of deviation 
due to secondary flow, i.e., 

AaAi* ( A <3 ) 5 Q - ( <5 - 6 eq) - { <$ “ 6eq)r. t g, where (.<5 - fieqjgg 

♦Numbers in brackets denote references at the end of the report.' • 
'Symbol definition in Appendix A. 
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is the midspan value. This has been clone; the results for each 
stator combined with various rotors are given in Appendix B. 

Figure 2 shows the difference between rotor and stator 
span wise variation of deviation. There the spanwi.se variation 
due to secondary flow is plotted for rotor 2 3B/ stator 20 
operating together [1] using the method described above. It can 
be seen that the rotor has overturning at the hub and 
■underturning at the tip with little variation between. This is 
probably due to the effect of rotation on the fluid in the blade 
wake. This fluid is centrifuged to the tip region. The gap at 
the tip would further increase the under turning due to the 
accumulation of low energy fluid at the tip. 

For stators with no end-wall gaps, there are two distinct 

K 

I maxirnurns of underturning that occurs near the end-walls as seen 

'• -in' Figure 2'." These maxirnurns are caused by the "rolling uo"-.of 
the two major secondary, viscous corner vortices that form 

I because of the passage pressure gradient. The wall side of the 

! vortex enhances overturning and the mid-passage side underturning 

ft 

j due to the direction of rotation, see Figure 3. 

1 

1 The physical parameters that can affect stator spanwise 

I 

•j distribution of deviation are blade camber, <j), blade aspect 

ratio. Ah, blade channel aspect ratio , AR c i solidity, • <r> 
end-wall boundary layer thickness at the inlet, represented by 
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the inlet boundary layer displacement thickness, d J , and 
perhaps hub/tip radius ratio. All of the above parameters 
are available for all the rotor/stator combinations except 
displacement thickness. However, Reference 3 gives values of 
o l for seven (7) of the rotor/stator pairs. A list of these 
combinations, along with pertinent physical parameters, is given 
in Table 1. 

A perusal of the data plots in Appendix B shows that, for 
most stators, there is a maximum underturning that occurs within 

. ' 10-20% of the end-walls With an unsy mmetr ica 1 exponential type 

variation between there and the end-wall or mid-span. Therefore, * 

r 

the variation .of ••• secondary flow induced deviation has been 
modeled 1 n four ( 4 ) parts: ( 1 ) the magnitude of • the 

maximum underturning, ( 2 ) the location of the maximum 
underturning, (3) the value of underturning or overturning at. the 
end-wall and (-4 ) .the variation between. A schematic of this 
model is shown in Figure 4. 

• . THE MAGNITUDE OF MAXIMUM UNDERTURNING 

c. . 

A series of physical models for maximum underturning that 
reflect various combinations of the parameters listed in Table 3. 
is given . in Table. 2. The criterion adopted for an acceptable 

kp 

||5 
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model is the ability to correlate A max within +1°. This value 
was chosen because that is the expected accuracy of the data [1.]. 
Example plots of data from Table 1 using some of these models are 
given in Reference 4. Hub/t.ip radius ratio was not included in 
these models because it was constant at 0.0. 

Model 1 was an attempt to account for the maximum 
underturning using commonly-avai lable geometrical parameters. 
The plot of Model 1 versus A max indicates that something 
essential has been left out, see Figure 5. A perusal of the 
literature on secondary flow [5, 6] indicates that the end-wall 
boundary layers have a significant effect on secondary loss and 
turning in cascades and rotors. This is natural, since secondary 
flow is the redistribution of the end-wall boundary layer as it 
passes through the blade row. The stators of Table 1 are those 
for which the boundary layer properties were available at the 
•tip.' • -Inlet -displacement' thickness, was chosen as the parameter .'to 
represent the b o u n d a r y layer e f feet. Figure 6 shows t h e 
variation of A with <5?. As can be seen, the correlation is 
significant . 

T h o r emainder of t h e m o d e 1 s t e s t e d h a d displace m e n t 
thickness in the numerator. For the thirteen remaining models, 
only 5 and 7 correlated all the data within +1°. Of these two, 
Model 5 has the best fit with the data. Figure 7 shows the 
linearized version of Model 5. From this, the relation between 


\ 
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the maximum value of underturning and the important physical 
parameters can be deduced: 


r 0 ( 6 *) 2 -. 3 / 4 

A max * — =— (deg) (1) 

L AR C * o J 


Equation 1 is a specific relationship modeling NASA middle 
'stage data [1]. However, the model should be applicable to 
similar subsonic compressor stages. Therefore, it would be 

useful to have Equation 1 in a general form. This can be done by 
•appropriately normalising the parameters in the bracket. The 
terms in the denominator are already dimensionless, so the camber' 
and displacement thickness are the parameters that should be 
modified. The camber can be non-dimensional i 2 ed by specifying 
its value in radians. It is physically logical to normalize the- 
displacement thickness by the mean blade height. If an inlet 
boundary layer is large compared to the blade height, the effect 
on spanwise deviation should be large and vice versa. A drawback 
in this case is that the data used to make up the model all had_ 
the same blade height of 10.16 cm and there is no acceptable 
data available with a different height. Therefore, the general 
model will be unproven for different blade heights. 



The non-dimensional linear model for maximum underturning is 
given below and in Figure 8 with the constant of proportionality 
adjusted to reflect camber, 0 , in radians and inlet 
displacement thickness normalized by blade height or span, 5^/s: 


- r0(<5j/s) 2 -| 3;/4 

Arr ... = 2.5 x 10 4 (deg) (2) 

" I. AR C - cr J 


LOCATION OF MAXIMUM UNDERTURNING 

Table 3 lists the location of maximum underturning for 
twelve (12) NASA rotor/stator combinations and the exit of Stator 
3 of the General Electric Low-Speed Research Compressor (7). A 
plot of this data, shown in Figure 9, shows that the maximum 
under turning . location, clusters between ten to twenty percent'. of 
span (s m '= 0.1-0. 2) from the hub or tip. A careful study of . the 
data indicates that there is no combination of parameters that 
model this result. All the NASA data is from test staged with a 
hub/ tip radius ratio of 0.8, so this could be a correlating 
factor, although the GE stator has a radius ratio of 0.7 and' it 
clusters with the NASA data. When the final model was compared 
to t h e d a t a i t w a » found that the b e s t f i t o c c u r red 
when the location of maximum underturning A rnax 7 was taken -as 
12.5 percent of the blade span (i.e., height) from the end-wall. 
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This agrees well with the average of the locations shown in 
Figure 9. Therefore, until additional data is aval 1 abl e , * the 
location of maximum underturning will be taken between ten and 
twenty percent of span from the hub or tip with an arbitrary 
value of twelve and one half percent, i.e.: 

®m-h = s m -.t = 0.125 (3) 


END WALL VALDES 

To completely model the spanwise variation of deviation, the 
value of under or overturning, A w , must be estimated at the 
wall. This can be done by extrapolating the data of Appendix B. 
The method used was to fit a curve through the 15%, 10% and 5% 
span data points and on to the wall. This is illustrated in 
Figure 10, where it can be seen that there is a region of 
"reasonable" maximum to minimum extrapolation. 

Since the value of deviation difference at the wall is the 
result of the same corner vortex that causes maximum 
u nderturning, and t h e strength of this vortex will be 
proportional to the difference between A max an<3 the value at 
the wall, the same parameters that modeled A max should model 
the difference: 


8 




_ _ 0( <5 */s) 2 

( ^max“ ~ 

AR c* tr 


Table 4 shows the appropriate parameters, with values of A w 
adjusted within the region of maximum to minimum extrapolation 
to best fit the model. Figure 11 shows that with the exception 
of one data point there is a linear relation between deviation 
difference and the vortex parameter: 

_ _ rlj(6p's) 2 ] 

A max “ - 5.7 x 10 5 (deg) (4) 

L AR C - 0- J 

Once A, nax is known, Ay can be calculated using Equation 4. 

SPANWISE VARIATION 

After the values of the m a x i m u m and w a 1 1 deviation 
difference and the location of the former have been found,. the 
variation between the mid-span and the wall must be determined. 
The variations plotted in Appendix A .indicate that a modified 
normal distribution, could be used. After a data fitting process 
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the following relations have been adopted: 

1) The variation between the location of maximum 
underturning to the mid-span (where A = 0) is 

modeled by: 

2 

A * A n , ax (e- (Gx/£) 1 (deg) < 5 > 

where x is the spanwise distance from the 
location of 5 max to the mid-span and s is the 
mean blade span or height,, 0 < x < 0.375 . 


The variation between the location of maximum 
underturning to the wall is modeled by: 


2 

a X i A - A ) [e"( 20 x/s) ] (deg) (6) 

A =• A max " 1 A max w' 1 


where x is the spanwise distance from t he 

location of A max to the wall and s is 

the mean blade span or height, 0 ~ 0.125. 

s 

Equations 5 and 6 apply in both the hub or tip region. 
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RESULTS AND DISCUSSION 

The model, Equations 2 through 6, has been compared with 
experimental half-span tip region distributions for the stators 
of Table 1 (using Appendix B). The results are shown in Figure 
12. There it can be seen that the prediction is fair to good for 
all stators with the exc f cion of stator 22 (Fig. 12f). 
Therefore, the present model can be considered to be a good 
approximation for subsonic NASA middle stages of medium hub, tip 
radius ratio. 


Before the model is applied outside the range of. the data 
base .[1], it should be validated against appropriate data, i.e., 
that which includes the inlet boundary layer displacement 
thickness. Unfortunately, additional data of this type is not 
currently available to the authors. Until such data, is 
available, the model should be treated as preliminary. 
However, there is one factor that has been demonstrated: the 
effect of the end-wall boundary layer on the spanwise variation 
of deviation. It has been long known that the end-wall boundary 
layer is a primary factor in the development of secondary flows.. 
However, most attempts to model the loss and deviation variation 
caused by end-wall effects have ignored the specification, 
estimation, or calculation of end-wall boundary layer thickness 
variation through the compressor. This can no longer be’ 



m 

fe * 


neqlocted. Furthermore, future research compressor data. should 
specify hub and tip boundary layer development for at least 
design Point operation. Such data would be of great benefit to 
the model analysts and computational fluid dynamicists. 

CONCLUSION 

It has been demonstrated that the following five (5) 
equation models can predict the spanwiso variation of deviation 
to a good approximation for NASA .subsonic . middle stages (applies 
for the hub and tip regions): 


Va 1 u o o f ma x i mu m 
under turning , 


,4rna-x. 


2 . 5 x 1 0 


-.4 


0( 6 j/s) 2, 3 


• AR C '* o' 


3/4 


(deg) 


Location of maximum under turni ng , 


(s) 


ma x 


0. IPS 


1 




Value at the wall, A^ * A, nax - 5.7 x 10 s 


Q( 6 f/s) : ’-l 


AR C * a - 


12 


Variation from location 
of maximum underturning 

to midspan# _ „ 2 

A = A m£ix ( 0 ~ ( 6x/u) ] (deg) 


Variation from location of 

maximum underturning __ 2 

to the wall. A* A, nax - ( A max - A w ) (e “( 20 J (deg) 


This model has not been tested outside the data base. However, 
it should prove useful for making engineering predictions of 
span wise variation of deviation for stators if ( 1 ) the 
distribution of two-dimensional blade element deviation is 
available, (2) there is a region of axi-symmetr ic core flow and 
(3) the inlet Mach number is less than one. 
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APPENDIX A 
SYMBOLS 

AR ~ aspect ratio 

AR C - channel aspect ratio 

s - span or blade height 

s m - fraction of the span to the location of 

maximum underturning from hub or tip 
x - distance along the span 
® - deviation angle 

A - difference, in. deviation , ( (,- 6 eq) - ( <5 - <5 eq) 5Q 
<5* - boundary layer displacement thickness 
cr ~ blade solidity 

^ ■* blade camber 

Subscripts 


1 

i n 1 e t 

so - 

in id- span 

oq 

equiva lent 

It 

hu b 

max “ 

maximum 

t 

tip 


w 


end-wa 1 1 




APPENDIX B 


Variation of Deviation Due to Secondary Flow for NASA Middle 
Stage Stators Combined with Various Rotors 














































•co bssss « T2^4£n2a 


























































KCUrrsL « ES8£ 

































































| , a* ^pj^g wBfeMfiffiM8j3B& >> wi8»4PBBw M BBw8Bwiillwi^ 


*J% 


ORIGINAL PAG2i W 
OF POOR QUALITY 


TABLE 2. MODELS FOR MAXIMUM UNDERTURNING 


Model No. 


1 


2 


3 



5 


6 


7 


Model, A = 
’ max 


K£ 

AR*AR„ /a 
c 


K<j)<$ i 
AR /a 


K<J> 6 i 
ATT“AR~ 


K$<5j 

AH* AR C • a 


K' R'S *) ? ' 
™c' a 


K-HSO 2 

:(AR c .-.o)-y f . . . . . 


K <J>(6*) 2 ' 5 

AR ♦ a ~ 
C ' 


Model No. 


8 


9 


10 


11 

12 


13 


14 


Model, A 


max 


KjCjll 

/M^o 


Kjfeiii ) 2 

AR C /a 


K<() (6 i ) 

Aft 


K<K«i V 


•k 

K(j) <5 i 

“7CR™" 


K(Si) 


K 6, 

w c ~r 


Nom e ncla t ure : AR - Blade Aspect Ratio 

AR C - Qiannel Aspect Ratio 

K - Constant 
a 

5i - End-Wall Displacement Thickness at Inlet 


<p ■ - Blade Camber at hub or tip 



TABLE 3. 


S PANWISE POSITION OF- MAXIMUM UNDERTURNING FOR NASA STATORS 


Stator 

— t 
s ia, h 

s m,t 

* 

5 i, t 

(in. ) 

ar : 

c 

AR 

°t 

: 

_ _ * 

! Stator s u s* Si,t 

1 (in.) 

AR c 

AR 

°t 

20 

(R238) 

.125 

.100 

. 020 

1 .90.' 

1.04 

1.16 


210 in u 

(R25D) ’- 1 

! 

1.24 

1.24 

0.9 

20 

(R24A) 

.100 

.150 ; 

. 027 

1.90; 

1.04 

1.6 


(R27D) * u:> 

; 

! 

1.24 

1.24 

0.9 

20C 

(R23D) 

.300 

.150 

. 018 

2.34-. 

•1.04 

2.0 


22 15 

(R2SB') 

1.63 

0.82 

1.8 

20C 

(R2.4B) 

.150 

.125 

. 032 

2.34- 

1.04 

2.0 


(R28D) * 15 ,2 ° * 01 

1.63 

. 0.82 

1.8 

2QB 

(R2SA) 

.150 

.100 

. 022 

1.42 : 

1.04 

1.2 ! 

: 


LSRC 

No. 3 .20 .10 

Stator 



do 

-if 70 

21 

(R26B) 

.150 

.100 

- 

2.48 

1.24 

1.8 

| 

Ay S; ? m,h = * 144 



T» O 
O -A 

O 

Si f 
O • 

c : 
~=-. . 

(R27A) 

.150 

.100 

— 

2.48 

1.24 

1.8 


Av s- s m>t - .121 



. < * 

21 

CR27C) 

..100 

. 150 

.027 

2.48 

.1.24 

KB 

i 

! 
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TABLE' 4 . TIP END-WALL VALUES OF THE 
DIFFERENCE IN DEVIATION 


r"" 1 " - 

■ 

1 

f Stator 

* 

"max j-t 
(deg) 

7 .. 

W, t 

u nax, t 

’ A w,t 

(deg) 

Symbol 


| 

! ' 20 
! (Rotor 23B) 

t ' - 

H 

°- 

■ • 

. 4,1 

■ 

O 

i 

i 20 

! (Rotor 2 4 A) 
I 

6,.s: 

• 

. 

■ 

-2.5- 

9.0 

. 

A 

| 

I 205 

j (Rotor 25A) 

- 

o. 6 

* 

-3.0 

■ 

% 6 

■ - - - 

□ 

i 

| 20C 

| (Rotor 23D) 
! 

2.4 

0.5 

1.9 

o 

! 

2GC 

(Rotor 248) 

m 

■■ 

6.6 



21 

(Rotor 27 C) 

6,0 

•1.5 

i mi 

( 

3 

22 

(Rotar 28D) 

1.2 

-4.8 

’ 

6.0 

• 

a 
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S PANWISE DEVIATION PERFORMANCE OF 1 

NASA STATORS IN COMBINATION 
WITH VARIOUS ROTORS AT MAXIMUM 
EFFICIENCY OPERATION 
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STATOR SPANW1SE VARIATION OF 
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FIGURE 4. • 

. MODEL FOR THE VARIATION OF 
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Nomenclature: A" - AS-(A<$)so» AS ■ 6-6eq 

A “Maximum Underturning 

max ’ 

(x\ . Fraction of span to maximum 

m ~ \s/ ' underturning from hub or tip 

h - hub 

t - tip 




FIGURE 5. 
NASA STATORS 

MODEL 1 
(deg) 





ORIGINAL PAUL UJ 

OF POOR QUALi »V 


FIGURE 6. 

MAXIMUM UNDERTURNING VERSUS 
INLET TIP DISPLACE M ENT THICKNESS 
FOR NASA STATORS 
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FIGURE 10. ILLUSTRATION OF REASONABLE REGION 

OF MAXIMUM TO MINIMUM EXTRAPOLATION 
FOR END-WALL VALUE OF DEVIATION 
DIFFERENCE 
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RMME - Region of Maximum to Minimum 
Extrapolation 
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